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Abstract 

We present the effect of Zn (S = 0) and Cu (S = 1/2) substitution at the Ni site of S = 1 Haldane 
chain compound Y2BaNiOs. 89 Y NMR allows us to measure the local magnetic susceptibility at 
different distances from the defects. The 89 Y NMR spectrum consists of one central peak and several 
less intense satellite peaks. The central peak represents the chain sites far from the defect. Its shift 
measures the uniform susceptibility, which displays a Haldane gap A ~ 100 K and it corresponds 
to an AF coupling J ~ 260 K between the near-neighbor Ni spins. Zn or Cu substitution does not 
affect the Haldane gap. The satellites, which are evenly distributed on the two sides of the central 
peak, probe the antiferromagnetic staggered magnetization near the substituted site. The spatial 
variation of the induced magnetization is found to decay exponentially from the impurity for both 
Zn and Cu for T > 50 K. Its extension is found identical for both impurities and corresponds 
accurately to the correlation length £(T) determined by Monte Carlo (QMC) simulations for the 
pure compound. In the case of non-magnetic Zn, the temperature dependence of the induced 
magnetization is consistent with a Curie law with an "effective" spin S = 0.4 on each side of Zn. 
This staggered effect is quantitatively well accounted for in all the explored range by Quantum 
Monte Carlo computations of the spinless-defect-induced magnetism. In the case of magnetic Cu, 
the similarity of the induced magnetism to the Zn case implies a weak coupling of the Cu spin to the 
nearest- neighbor Ni spins. The slight reduction of about 20 — 30% of the induced polarization with 
respect to Zn is reproduced by QMC computations by considering an antiferromagnetic coupling of 
strength J' = 0.1 — 0.2 J between the S = 1/2 Cu-spin and nearest-neighbor Ni-spin. Macroscopic 
susceptibility measurements confirm these results as they display a clear Curie contribution due 
to the impurities nearly proportional to their concentration. This contribution is indeed close to 
that of two spin-half for Zn substitution. The Curie contribution is smaller in the Cu case, which 
confirms that the coupling between Cu and near-neighbor Ni is antiferromagnetic. 

PACS numbers: 75.30.Hx, 75.40.Cx, 76.60.-k, 76.60.Cq 



I. INTRODUCTION 

One and two-dimensional Heisenberg an- 
tiferromagnetic (ID and 2D HAF) spin sys- 
tems are presently attracting strong inter- 
est due to the variety of ground states that 
these systems exhibit : spin-gap, spin-Peierls, 



RVB, superconductivity, etc. Effects of 
substitution in the HAF half-integer— and 
integer-spin chains,— ladders,— and planes^ 
have been studied both experimentally and 
theoretically to understand such ground 
states and their elementary excitations. 
Among these systems, one-dimensional 
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S = 1 HAF chains have attracted immense 
interest after Haldane^ conjectured that in 
case of integer-spin chains, there is a finite 
energy gap between the ground state and 
first excited state. The ground state of this 
system can be well described by Valence 
Bond Solid (VBS) model£ This model sug- 
gests that each S = 1 spin can be consid- 
ered to be a symmetric combination of two 
5=1/2 spins. With periodic boundary con- 
ditions, each S — 1/2 forms a singlet with 
the nearest-neighbor (n.n.) belonging to the 
adjacent site. Introducing a non-magnetic 
impurity, i.e. open boundary conditions, is 
equivalent to removing one S = 1. This is 
expected to give rise to two S = 1/2 at each 
end of the chain.— 

Y 2 BaNi0 5 (YBNO), an S = 1 one- 
dimensional HAF is a well established Hal- 
dane gap compound. YBNO can be de- 
scribed by the space group J mmm .— All the 
Ni 2+ ions are located at the center of highly 
compressed corner shared NiC>6 octahedra 
and form Ni - O - Ni chains along the a- 
axis (see Fig. 1 (a)). These octahedra are 
flattened along the chain axis and the Ni - O 
distances can be divided in two groups (Ni 
- 01 « 1.885 A and Ni - 02 « 2.186 A). 
Intra-chain Ni - Ni distance is equal to 3.77A 
(Ni 2+ ions interact via the oxygen ion), while 
the smallest inter-chain Ni - Ni distance is 
close to 6.6 A. Also, since these chains are 
separated by Ba and Y cations, the physi- 
cal properties are expected to be predomi- 
nantly one- dimensional (ID). The superex- 
change coupling J estimated from magnetic 
susceptibility is J ~ 285 K.— The magni- 
tude of the Haldane gap A rs 100 K has 
been obtained from magnetic susceptibility 
measurements.— i 10 ^ 1 Inelastic neutron scat- 
tering experiments (INS) also reveal the ex- 
istence of a gap A « 100 K. 912 between 
the ground state and first excited state and 
\J±/J\ < 5 x 10~ 4 , where J± and J are 
the inter-chain and intra-chain coupling con- 
stants, respectively— This proves the strong 
ID character of the Ni chains. 89 Y NMR 



experiments on polycrystalline YBNO, 14 also 
show the presence of a gap A = 80 ± 5 K and 
an anisotropic Knight shift which is respon- 
sible for the asymmetry in the NMR spectra. 

In-chain Zn 2+ , Mg 2+ (S = 0) substitutions 
at the Ni 2+ (S = 1) site of YBNQ W 6 ^ 
have been attempted. Bulk susceptibil- 
ity measurements confirm that non-magnetic 
substitutions break the chains and give rise to 
chain-end free spins.— Low-temperature spe- 
cific heat measurements 1 ^ on different doped 
samples have been quantitatively explained 
by an effective model,— which describes the 
low energy spectrum of non-interacting open 
S — 1 chain. Tedoldi et al^ carried out 
89 Y NMR experiments in 5% and 10% Mg 
doped YBNO for temperatures T > A. They 
observed multiple-peak 89 Y NMR spectra at 
different T and determined the T-dependence 
of the staggered magnetization near the im- 
purity from the shift of these peaks at dif- 
ferent T. From the site-dependence of the 
staggered magnetization, the characteristic 
length scale of its decay was found to be sim- 
ilar to the correlation length £(T) of the pure 
S = 1 HAF chain. At lower temperatures, 
the NMR spectra were found to broaden sig- 
nificantly and hence various features could 
not be resolved. These results prove that the 
non-magnetic impurity- induced effects reveal 
the intrinsic characteristics of the correlation 
functions of these systems, both their expo- 
nential shape and their extension. Such re- 
sults were quantitatively explained by Monte 
Carl o 19i20 | i2:L and Density Matrix Renormal- 
ization Group simulations . 21 > 22 i 23 

In this paper, we propose to extend 
such local studies to other impurities, both 
non-magnetic and magnetic and to greatly 
improve the quantitative aspects of these 
results. We, therefore, present a study 
of the susceptibility and 89 Y NMR spec- 
tra on Y 2 BaNii_ x Zn x 05 (YBNO:Zn) and 
YsBaNi^CuzOs (YBNO:Cu) (0 < x < 
0.02). Dilute concentration of substitutions 
were preferred in an attempt to observe well 
resolved impurity- induced effects on the 89 Y 
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NMR spectra far from the impurities and 
in a temperature regime T < A. We de- 
tail the experimental conditions in Sec. II 
and NMR and bulk susceptibility measure- 
ments on YBNO:Zn and YBNO:Cu will be 
presented in Sec. III. Zn and Cu substitu- 
tions are found to induce low-T Curie-terms 
in the bulk susceptibility, while 89 Y NMR 
spectra show a multi-peak structure. The 
multiple peaks are associated to an impurity- 
induced short-range staggered magnetization 
which appears very similar for both Zn and 
Cu. As discussed in Sec. IV, our results on 
YBNO:Zn prove the universality of the ef- 
fects of spinless (Mg and Zn) impurities in the 
S = 1 chain. These results are quantitatively 
reproduced by our QMC simulations. Effects 
of an S = 1/2 impurity (Cu) are contrasted 
with those of Zn substitution at various tem- 
peratures. Here again, QMC simulations fit 
the data well and enable us to deduce the 
magnitude of the exchange coupling between 
Cu and n.n. Ni, which is found weak and 
antiferromagnetic. Finally, the contribution 
of impurity induced magnetization to the to- 
tal susceptibility is extracted from NMR for 
the first time and has been compared to the 
macroscopic (SQUID) measurements taken 
on the same set of samples. All results can be 
well explained in a common VBS-like frame- 
work both for Zn and Cu. 

II. EXPERIMENTAL DETAILS 

Polycrystalline samples of 

Y 2 BaNi 1 _ a ,Zn E 05 and Y 2 BaNi 1 _ a; Cu a ;05 
(0 < x < 0.02) were prepared by standard 
solid state reaction techniques. Stoichiomet- 
ric mixtures of pre-dried Y 2 3 , BaC0 3 , NiO, 
ZnO, and CuO (> 99.9% pure) were fired 
in air at 1200°C. Single phase samples were 
obtained after several rounds of grinding, 
pelletization, and firing. In an attempt 
to reduce the parasitic Curie term in the 
susceptibility, as discussed in Sec. IIIA, we 
fired the samples in vacuum at 400°C, in 
order to reduce the excess oxygen present 



in the sample. X-ray powder diffraction 
was carried out using a PHILIPS PW 
1729 powder diffractometer, using Si as an 
internal standard. A Cu target was used 
in the diffractometer with A a „ = 1.54182 
A. Lattice parameters were determined 
from the diffraction pattern using a least 
square fit method. The lattice parameters 
for the" pure" YBNO sample were found to 
be 3.763(1) A, 5.762(1) A, and 11.334(3) A 
along a-, b-, and c- axes respectively, which 
are in agreement with the earlier results. 8 
The lattice parameters remain almost un- 
changed after Zn and Cu substitution. The 
nearly identical effective ionic radii^i of Ni 2+ , 
Zn 2+ and Cu 2+ , together with the dilute 
level of doping may be the cause of the 
unchanged lattice parameters. 

Magnetization M measurements were per- 
formed function of temperature T (1.8 
K < T < 300 K) and an applied field H (0 

< H < 50 kGauss) using a Superconducting 
Quantum Interference Device (SQUID) mag- 
netometer. 

NMR measurements were performed by 
standard pulsed NMR techniques in an ap- 
plied field H = 70 kGauss at different tem- 
peratures (50 K < T < 325 K). The typ- 
ical width of the applied 7r/2 pulse was 
tw — 6 [is. The line shape was obtained 
by Fourier Transform (FT) of the NMR echo 
signals. The 89 Y NMR Knight shift K(T) = 
[v(T) — u re f] I v re j was measured with respect 
to a standard YCI3 solution iy{T) is the res- 
onance frequency of Y 2 BaNiOs at different 
temperatures and z/ re j = 14594.15 kHz is 
the resonance frequency of YCI3 in our ex- 
periment). At higher temperatures (120 K 

< T < 300 K), the spectral width of the ?r/2 
pulse was sufficient to reliably obtain the full 
spectrum. However, at lower temperatures 
(T < 120 K), where the multiple resonance 
lines are spaced far apart in frequency, we 
varied the irradiation frequency to obtain the 
full spectrum as the spectral width of the sin- 
gle Ti / 2 pulse was insufficient. 

The Quantum Monte Carlo simulations 
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were performed using the Stochastic Series 
Expansion method , 25 > 26 with generalized di- 
rected loop techniques . 27 ^ 28 ' 29 Calculations of 
magnetization profiles were done for long 
runs of about 2.10 8 Monte Carlo steps for an 
isotropic S = 1 Heisenberg model, with an 
imposed magnetic field of about 70 kGauss, 
corresponding to the field used in the NMR 
experiments. 

III. RESULTS 

A. Bulk Susceptibility 

Susceptibility xCO °f Y2BaNix_3.Zn3.O5 
and Y 2 BaNii__Cu_0 5 (0 < x < 0.02) was 
measured clS db function of T in an applied 
field of H = 5 kGauss (Fig. 2). We checked 
for the possible presence of spurious ferro- 
magnetic impurities by measuring M vs H 
isotherms at different temperatures in or- 
der to separate the linear (paramagnetic) 
and non-linear (ferromagnetic) contributions. 
The ferromagnetic contribution corresponds 
to about 2 ppm of ferromagnetic Fe impu- 
rity which amounts to at most 0.2% of the 
susceptibility. It is, hence, neglected in our 
analysis. The paramagnetic susceptibility for 
the "pure" sample is in rough agreement with 
the results published earlier The suscep- 
tibility decreases as e~ t J y/T with decreasing 
temperature, as a consequence of the Haldane 
gap.— At lower temperatures (T < 50 K), 
there is a Curie- Weiss like C/(T + 9) upturn 
in the susceptibility. This parasitic Curie 
term in the "pure" sample might presumably 
result from various sources such as natural 
chain breaks, presence of Ni 3+ ( result of 
small oxygen non-stoichiometry) , and extrin- 
sic paramagnetic phases present in the sam- 
ple. A substantial reduction in this low-T 
contribution to the susceptibility (i.e. reduc- 
tion in the Curie term C) was achieved by 
vacuum annealing the powders at 400°C, as 
evident from the inset of Fig. 2 (a). A sim- 
ilar reduction in the parasitic Curie-term in 
a S = 1/2 chain compound, Sr 2 Cu03, was 



observed earlier by Ami et al. 32 after anneal- 
ing in inert atmosphere. Reduction in the 
Curie term in our samples might be due to 
removal of excess oxygen leading to a reduc- 
tion in the amount of Ni 3+ due to the above 
process. A discussion about intrinsic defects, 
effect of vacuum annealing and the implica- 
tions on SQUID and NMR measurements is 
carried out in Sec. IVD. All measurements 
with the substituted samples reported here- 
after are on the vacuum annealed samples. In 
case of the substituted samples, the suscep- 
tibility at lower temperatures increases with 
increasing dopant content (shown in Fig. 2 
(a,b), Inset), which indicates an increase in 
the Curie-term for the substituted samples. 

B. 89 Y NMR 

In YBNO, each 89 Y has two n.n. Ni 2+ 
from one chain and two next-nearest- 
neighbor (n.n.n.) Ni 2+ from two different 
chains (Fig. 1 (a, b)). However, 89 Y is ac- 
tually coupled to the Ni 2+ spins belonging to 
n.n.n. chains, not to the n.n. one. Indeed, 
there is no appreciable exchange between Y 
and Ni 3d or oxygen p orbitals belonging to 
n.n. chain. In contrast, there is an appre- 
ciable overlap between the n.n.n. -oxygen p n 
orbital with the 5s- orbital of Y 3+ and hence 
with the n.n.n. Ni - O - Ni chains. 

We performed 89 Y NMR on the "pure" 
and substituted samples at different temper- 
atures. In case of the annealed and non- 
annealed "pure" samples, we observed one 
asymmetric central peak (inset: Fig. 3) 
shifted from u re f in the Fourier transformed 
spectrum, together with less intense peaks 
(satellites) on both sides of the central peak 
(see Fig. 3). The shift in the 89 Y NMR 
peak of the "pure" YBNO arises from the 
coupling of Y nuclei to Ni 2+ spins of two 
different n.n.n. chains. Since our mea- 
surements are on randomly aligned polycrys- 
talline samples, asymmetric lineshape of the 
central peak (mainline) in the annealed sam- 
ple corresponds to a powder pattern due to an 
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asymmetric hyperfine coupling tensor ^U/ 14 
and anisotropic susceptibility. 11 The vacuum 
annealing procedure leads to a large narrow- 
ing of the mainline and a decrease of the 
satellites intensities. The amount of defects 
corresponding to these satellites can be de- 
termined by measuring their area. Using an 
asymmetric Lorenzian fit, the fractional area 
of one of these satellites decreases from 3.2% 
to 0.8% after annealing. These satellites in 
"pure" YBNO presumably originate from the 
in-chain defects, which will be discussed in 
details in Sec. IVD. 

In substituted compounds, the shift of 
the mainline for the substituted samples, i.e. 
Knight shift K(T), remains unchanged from 
the "pure" sample at different temperatures 
(Fig. 4) and compares well with the results 
published earlier.— This implies that the Hal- 
dane gap remains the same after impurity 
substitution. In contrast, the introduction 
of impurities (Zn/Cu) at the Ni-site of the 
chain, results in a change in the local envi- 
ronment of the impurity. This change is evi- 
dent by the appearance of satellites on both 
sides of the mainline. For example, in Fig. 
5, we compare the NMR spectra of "pure", 
0.5%, and 2% Zn substituted samples at 120 
K (Fig. 5 (a)) and 0.5% and 1% Cu substi- 
tuted samples at 300 K (Fig. 5 (b)). The in- 
tensity of the impurity-induced satellites in- 
creases with increasing impurity concentra- 
tion keeping their shift from the mainline 
constant. This establishes that the satellite 
peaks are a direct consequence of substitu- 
tions and their positions are independent of 
their concentration , i.e. of the length of the 
chains. Furthermore, fractional intensities of 
these Cu and Zn-induced satellites are sig- 
nificantly larger than the weak satellites ob- 
served in the annealed "pure" sample. Hence, 
they can safely be attributed to Zn or Cu 
neighbors. Each of these satellites represents 
89 Y nuclei at a given distance from the im- 
purity. Since one Y nucleus is coupled to two 
different n.n.n. chains (see Fig. 1 (b)), for a 
concentration x of the impurity, there should 



be Ax Y-nuclei contributing to each satellite. 
Using an asymmetric Lorenzian fit (Fig. 6), 
we extracted the fractional areas of one satel- 
lite peak to the whole spectrum (outermost 
satellite on the right side of the mainline, 
marked by 1 in Fig. 6). For a nominal con- 
centration x = 2% of Zn or Cu, our fit leads 
to an effective concentration of (1.9 ± 0.1)% 
for Zn and (1.8 ± 0.2)% for Cu. This good 
agreement between the nominal stoichiome- 
try and the actual in-chain number of defects 
shows that the substitution amount is well 
controlled here. 

Comparison of Cu and Zn 89 Y NMR spec- 
tra done in Fig. 7 shows that Cu induces 
a staggered magnetization like Zn. But the 
difference in the shifts of the NMR satellites 
shows that the amplitude of the induced mag- 
netization near Cu is smaller than that for 
Zn. Representative 89 Y spectra are shown 
at different temperatures in Figs. 8 (a) and 
8 (b) for 2% Zn and 0.5% Cu substituted 
samples respectively. The Ni 2+ ions closest 
to the impurity are the most affected by the 
substituent and give rise to the most shifted 
satellite indexed as the I s * neighbor (1 in 
Figs. 8 (a) and 8 (b)). Further neighbors are 
numbered in accordance with the decreasing 
shifts of satellites with respect to mainline. 
Satellites with decreasing shifts alternate on 
each side of the mainline. This shows that 
the corresponding Ni sites have a magnetiza- 
tion, which alternates in sign in addition to 
the non-perturbed magnetization represented 
here by the mainline shift. The appearance 
of such induced staggered magnetization near 
the defects is indeed expected by theoretical 
computations , 21 ^ 3 as will be quantitatively 
detailed hereafter. 

In the case of Zn, the nearest Y nuclei to 
the impurity (marked by in Fig. 1 (b)) 
couple to both Zn itself and one undisturbed 
chain. Therefore, its corresponding NMR 
line should be shifted by half of that of the 
mainline with respect to v re fi as Zn is non- 
magnetic and hence does not contribute to 
K. Furthermore, this satellite should be half 
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the intensity of the other satellites, from its 
specific position. In fact, this satellite has 
been observed clearly in Fig. 8 (a) where 
it is labelled as "0": it is half as shifted as 
the mainline and half as intense as the other 
satellites. In the case of YBNO:Cu, such a 
peak is not observed as seen in Fig. 8 (b). 
This is not surprising as the 5=1/2 Cu-spin 
will now contribute to its shift, which should 
now become large, and would be severely af- 
fected by a large spin-lattice relaxation rate. 
Therefore, the absence of such a peak is a 
clear local evidence that Cu has indeed suc- 
cessfully substituted at a Ni site. 



The shifts of the observed satellites were 
measured for different temperatures and con- 
centrations. We note the shifts of the i th 
satellite as 5ui = Ui — v main , where Vi and 
Vmain are frequencies corresponding to its 
center of gravity and that of the mainline re- 
spectively. Due to the asymmetric powder 
pattern shapes, consideration of the centre of 
gravity for each peak yields a better descrip- 
tion of the impurity-induced effects in the 
chain. In Fig. 9 (a), we present these shifts 
of the satellites (5vi) at various temperatures 
for the spinless Zn impurity. On the same fig- 
ure, we plot as well the data published earlier 
for the Mg impurity, which is found identical 
to that of Zn. 17 Therefore, the induced mag- 
netization by both Zn and Mg is quantita- 
tively the same. In the present study, the use 
of samples with smaller impurity content and 
less natural defects allowed us to largely im- 
prove the resolution and to detect more satel- 
lites in a wider temperature range. A similar 
plot is shown in Fig. 9 (b) for YBNO:Cu. 
These measurements allow us to get a direct 
measure of the amplitude and spatial shape 
of the magnetization induced by both non- 
magnetic and magnetic impurities, as will be 
discussed later. 



IV. DISCUSSION 

A. Determination of the parameters 
of the Haldane Hamiltonian of the pure 
system 

The bulk susceptibility xCO i n the "pure" 
YBNO contains the contribution from some 
unknown defects and extrinsic paramagnetic 
phases in addition to the intrinsic gap con- 
tribution. On the other hand, in 89 Y NMR 
spectra, defects and intrinsic behavior can be 
easily distinguished since defects broaden the 
mainline or induce satellites, whereas the in- 
trinsic behavior is measured by the mainline 
shift. Hence, the T-dependence of the "gap 
term" is best reflected by the temperature 
dependence of the NMR Knight shift K(T). 

Indeed, the Knight shift 89 K(T) of the 
central line is expected to vary linearly with 
the uniform susceptibility Xu(T) through 

2 A 

K(T) = — $-x«(T) + K(0) (1) 
HbN a 

Here, is the 89 Y - Ni Hyperfine cou- 
pling constant and is independent of temper- 
ature, N a is the Avogadro number, and K(0) 
is the chemical shift independent of temper- 
ature. The NMR shift K(T) results from 
two n.n.n. chains, while Xu{T) represents the 
intrinsic susceptibility of one S = 1 chain, 
which explains the factor 2. As reported in 
Fig. 4, the Knight shift K(T) of the cen- 
tral peak decreases with decreasing T and 
remains unchanged for all the samples at the 
corresponding T. The observed exponential 
decrease in K(T) at low-T is direct evidence 
of the presence of the Haldane gap. Our 
measurements on different samples at differ- 
ent temperatures match well with earlier re- 
sults of Shimizu et al. for "pure" YBNO^ (as 
shown in Fig. 4). This indicates that the gap 
remains unchanged for impurity substituted 
samples. I. Affleck showed that the low tem- 
perature part of the susceptibility should fol- 
low a/ A/Te~i\— Previous studies used this 
model to fit the low-T part of both NMR 
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and macroscopic susceptibility measurements 
in order to extract the Haldane gap and the 
superexchange coupling J. This model also 
fits our data, as shown by a solid line in 
Fig. 4 for the CuO.5% compound, leading 
to K(0) = 340 ppm and A = (88 ± 7) K, 
which is in good agreement with the results 
published earlier . 13 > 14 

However, we propose here a more rigor- 
ous way to determine J. We used QMC sim- 
ulations to directly fit our data in the full 
temperature range. We have calculated the- 
oretically, by QMC simulations, the uniform 
susceptibility Xu versus TJ J i.e. the suscepti- 
bility of the infinite 5 = 1 chain theoretically, 
by QMC simulations. The Knight shift K ob- 
tained from our NMR experiments has then 
been plotted versus this \u for different J val- 
ues in Fig. 10 with T as an implicit parame- 
ter. The T — > value of K is constrained 
by Shimizu's NMR measurements,— which 
allowed to accurately determine m K(T — > 
0) = 370 ± 10 ppm, a value compatible with 
our measurements. Using this constraint, a 
good linear fit is obtained for J = 250 ± 10 
K (shown by the solid line) . This proves that 
the uniform susceptibility measured by NMR 
is well accounted by QMC simulations up to 
T = 350 K. This procedure of determining J 
allows for a much sharper constraint on the 
determination of J than the previous deter- 
minations. We believe that this procedure 
enables the determination of J more reliable 
than before. The corresponding theoretical 
Haldane gap^ would be A = 0.41 J = 102 ±4 
K, a value compatible with our first rough es- 
timate and with the previous experiments. 11 
This fitting method also allows a determina- 
tion of the hyperfine coupling Ahf indepen- 
dent of any macroscopic susceptibility mea- 
surement. We performed similar fitting for 
all substituted and pure YBNO shifts plotted 
in Fig. 10. This leads to A hf = (4.7 ± 0.2) 
kGauss//i£- This value should be more re- 
liable than earlier determinations done by 
Shimizu et a/.— (Ahf = 4.9 kGauss//xe) and 
by Tedoldi et a/.— (Ahf = 6.5 kGauss//^), 



which relied too much on the actual measure- 
ments of the susceptibility of the pure com- 
pound, as this measurement is always ham- 
pered by the large low-temperature Curie 
contributions. 

B. Impurity induced staggered mag- 
netization 

1. Spatial dependence of the induced stag- 
gered magnetization 

In Figs. 11 (a) and 11 (b), we depict the 
site-dependence of |£z/j(T)| at different tem- 
peratures for YBNO:Zn and YBNO:Cu, re- 
spectively, on a logarithmic scale. At a given 
T, the observed linear behavior shows that 
\8ui(T)\ varies exponentially function of 
distance from the impurity. In fact, the shifts 
can be fitted at all temperatures for both Cu 
and Zn by: 

\Sui (T)\ = \8v 1 (T)\e~^^ (2) 

Here £,i mp (T) is the extension of the impurity- 
induced staggered magnetization. The values 
of ^imp(T) extracted from our measurements 
are represented on Fig. 12 and compared 
with those obtained by Tedoldi et al. in the 
case of Mg substitution.— It clearly appears 
that Cu, Zn and Mg induce a staggered mag- 
netization with the same spatial extent. Fur- 
thermore, this spatial extent £j mp (T) is found 
almost identical to the correlation length 
computed for an infinite 5 = 1 chain with no 
defects.— The experimental values of £j mp (T) 
are systematically found about 5% higher 
than the values of £ (T)^ but this discrepancy 
is still within experimental accuracy (e.g. at 
200 K, f imp = 2.5 ± 0.2, while f = 2.3). The 
lowest temperature data, which were only 
taken in the case of Cu impurities, depart 
a little bit more from the value of £(T). The 
good agreement found shows that for both 
5 = and 5 = 1/2 impurities, the charac- 
teristic length of the staggered magnetization 
near defects reveals indeed the intrinsic £(T) 
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of the infinite S = 1 chain, in the explored 
temperature range. 

2. Comparison with QMC for Zn and Cu 

For both 5 = (Zn) and S = 1/2 (Cu) 
impurity, we compared our experimental re- 
sults to those predicted by QMC calculations 
for the local susceptibility around impuri- 
ties. As discussed earlier, the total magne- 
tization Mi at the i th site in the chain is the 
sum of the uniform magnetization and the 
impurity-induced magnetization i.e. Mi = 
gfiB < Si >= gfi B (< S z > u +5 < S l z >). 
We can determine directly, from the NMR 
data for 5vi and the mainline shift, the ratio 
of the impurity-induced staggered magneti- 
zation and uniform magnetization 5 < S\ > 
/ < S z > u , using the value of K(Q) = 340 
ppm deduced from Fig 4 and 10. This exper- 
imental quantity could be directly compared 
with that computed by QMC simulation for 
the spatial dependence of the spin polariza- 
tion at chain ends. 

In case of the S = impurity, we per- 
formed such QMC calculations considering 
J = 230 K, 250 K and 270 K at T = 
100 K. The computed impurity-induced stag- 
gered magnetization decays exponentially 
away from the impurity as in the experimen- 
tal results, following Eq. 2 (Fig. 13 (a)). The 
computed extension ^ imp matches well with 
the previous results^ considering J = 270 
K. Further, we found that the experimental 
ratio of the impurity-induced staggered mag- 
netization and uniform magnetization match 
well with QMC for J = 270 ±10 K as shown 
in Fig. 13 (a). This value of J is consis- 
tent with the one independently determined 
from the uniform susceptibility far from the 
defect J = 250 ± 10 K. This proves the co- 
herence of our analysis, as the same coupling 
J can account both for the susceptibility of 
the "pure" compound, and for the shape, am- 
plitude and extension of the induced mag- 
netization near spinless-defects. The deter- 
mination of Ahf enables us to convert the 



NMR shifts into absolute magnetization val- 
ues. Here, we have calculated < S\ >z n , 
from our experimental data using the relation 
< S\ >= 2tt(2z/; - (v ref + v main ))/2^gA hf . 
Here, 7/2tt = 2.0859 MHz/Tesla is the 89 Y 
gyromagnetic ratio and g is the Lande g- 
factor. The temperature dependence of < 
S l z > is also well fitted by the QMC simu- 
lations as exemplified in Fig. 13 (b). 

In case of the S = 1/2 substitution, theo- 
retical calculations were done earlier to pre- 
dict the behavior of S = 1/2 impurity in 
S = 1 chain.^iss The S = 1/2 impu- 
rity is modelled by considering that its cou- 
pling J' to its near-neighbor Ni can be dif- 
ferent from J. We have performed QMC 
calculations for various J' ranging from — J 
to +J. The magnetic impurity is found to 
induce a staggered magnetization as in the 
spinless case. Furthermore, its extension is 
found identical for temperatures higher than 
J'. These results are then in perfect agree- 
ment with our data. For T > J', the only 
changes, which begin to appear when com- 
paring the S = 1/2 and S = induced ef- 
fects in QMC results, is the amplitude of the 
staggered magnetization. We found that for 
a ferromagnetic coupling of the impurity to 
the n.n. Ni (J' < 0), the staggered mag- 
netization on the 1 st , 2 nd and further neigh- 
bors is larger than that for a non-magnetic 
impurity. In contrast, an antiferromagnetic 
coupling (J' > 0) leads to a smaller mag- 
netization. The experimentally observed Sui 
for YBNO:Cu is less than that for YBNO:Zn 
(J' = 0) at the corresponding temperature, 
as exemplified in Fig. 7. This strongly sug- 
gests that the Cu - Ni coupling is then an- 
tiferromagnetic. A quantitative comparison 
between QMC results and experiments en- 
ables us to deduce J'. In Fig. 14, we com- 
pared the experimental ratio {o~Vi)cu/{& v i)zn 
to that theoretically obtained for < 5S l z >c u 
/ < SS l z >zn for different values J' /J. The 
experimental and theoretical results match 
well for 0.1 J < J' < 0.2 J, which allows to 
deduce an estimate J' = (0.15 ± 0.05) J. 
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3. Determination of the total induced sus- C. Temperature dependence of the 
ceptibility macroscopic susceptibility 



The magnetization induced by impurities 
can be written as 8 Mi = g\iB& < S\ >, where 
5 < S l z >= 2ii5h'i/'ygAhf. This conversion en- 
ables us to deduce from NMR, the contribu- 
tion of substituted impurities xfrnp R to the 
macroscopic susceptibility. As the impurity 
induced magnetization decays exponentially 
function of site on either side of the impu- 
rity in the ID chain, Xi^ R can De expressed 
as: 



xi CO 



SMi(T) 



(3) 



X 



nmr (t) = 2xN AX i(T) ^ i _ 1 

imp ) ^qq /_^\ V 



(4) 

Here, Xi(T) is the impurity induced sus- 
ceptibility on the I s * neighbor to the impu- 
rity and Na is the Avogadro number, x is 
the amount of impurities per Ni in %. The 
variation of Xim^ R (T), shown in Fig. 15, can 
be fitted to x?m P R (T) = C NMR /T. From 
the fit, we get Cf n MR = (5.6 ± 0.6) x 10" 3 
cm 3 K/mole/% Zn and C^ R = (4.7 ± 
0.5) x 10" 3 cm 3 K/mole/% Cu. These Curie 
terms lead to the corresponding spin value S 
on each side of the impurity through C = 
2N A g 2 S{S + l)fj? B /3k B . For Zn, one finds 
an impurity-induced "effective spin" S = 
0.4 ± 0.06 on either side of the impurity in 
the chain, which is close to the VBS pre- 
dicted S = 0.5. In case of YBNO:Cu, as 
Xc%f * s a l so f° un( i to follow a Curie behav- 
ior, the impurity-induced susceptibility on ei- 
ther side of the impurity behaves almost as a 
free spin in the explored T-range. The lower 
value of C%™ R than Cf n MR is well accounted 
for, hereagain, by QMC simulations, if the 
Cu-spin couples weakly and antiferromagnet- 
ically to the nearest-neighbor Ni-spin. 



As already stressed, the macroscopic sus- 
ceptibility presented in Fig. 2 results from 
distinct contributions: the uniform contribu- 
tion Xu{T) which exhibits the Haldane gap 
which decreases exponentially with decreas- 
ing temperature; the intrinsic defects such as 
Zn or Cu, and the other unknown defects 
and paramagnetic species, which are already 
present in the pure compound. The latter 
two contributions mainly give rise to a Curie- 
like magnetization in the low-temperature re- 
gion. We plotted these different contribu- 
tions in the case of pure YBNO in Fig. 16. 
Note that the total susceptibility not only 
consists of the Haldane term and a Curie 
term, but also includes a residual non-Curie 
term, which could be due to some defects 
that are not purely paramagnetic. This small 
residual non-Curie term exists in all our sam- 
ples, as evidenced in Fig. 17 and does not 
show any obvious connection with the sub- 
stitution. It, therefore, limits any quanti- 
tative analysis of the substitution effect on 
macroscopic susceptibility. Furthermore, the 
large contribution of the Curie term present 
at low-T even in the pure compound also 
severely limits any analysis, contrary to NMR 
experiments for which the Zn and Cu effects 
could be isolated. However, at low temper- 
atures, the total susceptibility can be tenta- 
tively modelled by: 



(jSQUID 

Xsquid(T) = Xo+ T + e + Xu(T) (5) 

Here, Xo * s the constant term, which orig- 
inates from Van-Vleck paramagnetism and 
core diamagnetism, C Q UID is the Curie con- 
stant. As discussed in the previous sections, 
the determination of the superexchange cou- 
pling J = 260 ± 10 K allows us to determine 
independently Xu{T) by using the QMC sim- 
ulation results, which fit perfectly the main- 
line NMR data. In order to estimate the 



9 



Curie constant C and 9 reliably, we there- 
fore fitted Xsquid(T) for the substituted as 
well as the pure samples to Eq. 5 in the tem- 
perature range 2 K < T < 25 K. We then fit 
our susceptibility data (T < 150 K) for all 
the samples to extract the values of x keep- 
ing the value of C and 9 obtained from the 
low-temperature fit. From the fit (shown by 
solid line in Fig. 2 (a)), we found that x re- 
mains almost unchanged for all the samples 
(x ~ 0.3 x 10 -3 cm 3 K/mole). A small 9 pa 
1 K fits the data reliably for all the samples. 

D. Comparison of NMR and macro- 
scopic susceptibility in pure YBNO: na- 
ture of intrinsic defects 

Comparison of NMR and SQUID mea- 
surements in pure non- annealed and annealed 
sample allows us to get a better understand- 
ing of the various types of defects present 
in the "pure" YBNO. The annealing pro- 
cedure leads to a reduction of the Curie- 
term C SQUID from 25.4 x 10~ 3 cm 3 K/mole 
to 13.7 x 10~ 3 cm 3 K/mole. Such reduction 
shows that annealing cures part of the de- 
fects in the compound. The corresponding 
89 Y NMR spectra of the non-annealed and 
annealed compounds are presented in Fig. 3. 
In the non-annealed sample, satellite peaks 
observed at frequencies smaller than the cen- 
tral peak, are of higher intensity than those 
noticed on the higher frequency side. Af- 
ter annealing, the intensity of the satellites 
reduced on both sides of the central peak 
and the intensity of the left side satellites al- 
most equals to those in the right side of the 
mainline. Hence, there are probably various 
kinds of defects: in-chain defects which give 
rise to satellites exactly like Zn would do, 
and other defects which lead to the broad 
low-frequency feature. There must be also 
magnetic entities at other non-chain sites, 
which do not affect the 89 Y NMR signal, since 
SQUID data show a large residual Curie term 
in the annealed compound. The in-chain de- 
fects could be chain breaks, due to the poly- 



crystalline nature of the sample, vacancies at 
Ni-site and presence of Ni 3+ in the chain. The 
other defects could be Ni 3+ in the intersti- 
tial places. Presumably, annealing in vacuum 
reduces the quantity of chain breaks (some 
Ni-ion might migrate to the vacant sites on 
annealing at 400° C) and Ni 3+ in the sam- 
ple, which is responsible for the reduction 
of satellite intensity. Annealing results in a 
reduction from 0.8% to 0.2%/Ni of in-chain 
defects. If one assumes that each of these 
in-chain defects is giving rise, as in the Zn 
case, to two S — 1/2 spins at low-T, the 
reduction of 0.6% percent of defects in the 
chain by annealing should lead to a decrease 
in the Curie constant of 4.5 x 10~ 3 cm 3 K/mole 
(0.6 x 2 x 3.75 x 10- 3 cm 3 K/mole). Our ob- 
served reduction in Curie constant was higher 
than this by 7.2 x 10~ 3 cm 3 K/mole. So the 
annealing not only reduces the chain defects 
but might also reduce the number of Ni 3+ 
(present in YBNO) or remove some of the 
excess oxygen present in the other extrinsic 
paramagnetic phases (say, BaNi02+:r). Some 
of these defects are however, still present in 
smaller concentration, even after annealing. 

E. Comparison of NMR and macro- 
scopic susceptibility in YBNO:Zn and 
YBNO:Cu: the low-T regime 

A plot of 1/xsquid versus temperature 
done in the insets of Fig. 2, demonstrates 
that low temperature susceptibility is dom- 
inated by a Curie term, which increases in 
the substituted compounds. In Fig. 18, 
we report the corresponding Curie constant 

C SQUID with 

increasing impurity content for 
Cu and Zn. This proves that the substitution 
induces a paramagnetic contribution, as al- 
ready has been shown by NMR at higher tem- 
peratures. This contribution seems to persist 
at temperatures as low as T = 2 K. How- 
ever, the "pure" Curie term appears simi- 
lar in magnitude to the Zn or Cu contribu- 
tions for concentrations of about 0.5%. A lin- 
ear behavior of C SQUID with increasing Zn 
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or Cu concentration is expected. However, 
we found that the variation of C S ® UID with 
substituent concentration is non-linear and 
is steeper for higher concentrations. This 
limits any straightforward quantitative in- 
terpretation of the data. This limitation is 
mainly due to the fact that in substituted 
samples, a sizable part of C Q originates 
from defects other than the substitution itself 
(i.e. background Curie term as discussed in 
Sec. IVD). The precise variation of this back- 
ground Curie term with the substituent con- 
centration is unknown. It appears that im- 
purity substitution (in small amounts) at Ni- 
site helps to decrease this background Curie- 
term as C SQUID for YsBaNio.ggsCuo.oosOs is 
almost identical to that of the "pure" whereas 
NMR showed that Cu induces a paramag- 
netic contribution. A rough estimate of the 
impurity contribution can still be made us- 
ing the increase of C ® UID for high enough 
dopant content. A linear fit (shown by solid 
lines in Fig. 18) without taking into account 
the "pure" composition leads to Cz® UID = 
(6.4 ± 0.7) x l(r 3 cm 3 K/mole/% Zn and 
c sQuid = ^ 3 ± Q g) x 10 -3 cm 3 K/mole/% 

Cu. 

In the Zn case, the NMR measurements 
allowed us to extract C^' IR = (5.6 ± 0.6) x 
10~ 3 cm 3 K/mole/% Zn. There is a good 
agreement between macroscopic and local 
measurements even though they apply in dif- 
ferent temperature ranges (T < 30 K and 80 
K< T < 350 K respectively). This value is 
not that far from the VBS simplified picture 
where Cvbs — 7.5 x 10~ 3 cm 3 K/mole& origi- 
nates from the two chain-end 5=1/2 spins 
released by the introduction of Zn. 

In the Cu case, if the doped 5 = 1/2 
Cu acts clS db free spin, its susceptibility 
should add to that of the induced staggered 
magnetization measured by NMR, leading to 

x SQUID + x NMR {n AnyhoW; we 

found Xcu UID t° be smaller than Xcu IR = 
(4 - 7 ^ a5) x 10" 3 cm 3 /mole/% Cu. This indi- 
cates that the introduced 5=1/2 does not 
act as a free spin, but magnetizes in the an- 



tiparallel orientation to the induced chain- 
end moments, at least at low-T. This result 
is indeed expected if the Cu interacts anti- 
ferromagnetically with the n.n. chain-end Ni 
spins, in which case xcu < is expected for 
T < J'. Such conclusion confirms then the 
comparison done previously between QMC 
simulations and our NMR results, which also 
leads to an antiferromagnetic coupling be- 
tween Ni and Cu, with J' « 40 K. 



V. CONCLUSION 

We have probed impurity induced mag- 
netic effects in the Haldane chain compound 
YBNO via spinless (Zn) and 5=1/2 (Cu) 
substitutions at the Ni site. The Haldane 
gap remains unchanged upon substitution. 
The intensity analysis of the NMR satellites 
proves that the amount of dopant impurity 
present in the sample agrees well with the 
nominally intended stoichiometry. 89 Y NMR 
measurements show that Zn or Cu induce 
a staggered magnetization in their vicinity. 
Its magnitude decays exponentially with a 
characteristic length equal to the correlation 
length £ of the infinite chain either for non 
magnetic or magnetic impurities. The fit of 
the NMR shift data on a large T range for 
the pure system with the results of QMC 
computations allows us to deduce a reliable 
value J = 250 ± 10 K. For the Zn impurity 
the spatial dependence and magnitude of the 
staggered susceptibility observed by NMR 
matches perfectly the QMC calculated values 
with J = 270 ± 10 K. For 70 K < T < 350 
K, the total impurity induced Curie suscepti- 
bility in the chain corresponds to an effective 
spin 5 = 0.4 ± 0.06 on either side of the im- 
purity for the Zn doped case, which is close 
to the theoretically predicted 5 = 1/2 exci- 
tation at T = 0. 

In this high-T range, the induced stag- 
gered magnetization is slightly smaller for Cu 
which can be explained if the 5 = 1/2 Cu 
spin is coupled antiferromagnetically with the 
near neighbor Ni. A fit with the QMC cal- 
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culated spin polarization allows us to deduce 
the magnitude J' = (0.15 ± 0.05)J of this 
AF coupling. This occurrence of an AF cou- 
pling between the Cu and the Ni is confirmed 
by the fact that the macroscopic susceptibil- 
ity induced by Cu as measured by SQUID is 
even smaller for Cu impurities than for Zn 
impurities for T < 30 K. This confirms that 
in this T-range the Cu magnetization opposes 
that of the local moments of the chain-ends, 
as can be expected for an AF J'. 

Further experiments are needed to follow 
the low-T behavior of the local magnetiza- 
tion both for Zn and Cu. In the case of Zn 
this would allow to check whether a Curie 
law applies to the lowest temperatures, be- 
low T « A, as suggested by macroscopic 
measurements, a matter still debated theo- 
retically. As for the Cu case this would al- 
low us to better follow the modifications of 
local magnetization for T « J' as strong 
deviations from that observed for spinless 
impurities are expected. In case of high-T c 
cuprates (e.g. YE^CuaOy), substitution of 
Cu by magnetic (Ni) or nonmagnetic (Zn) 
impurity— generates an induced polarization, 



with same temperature dependence of the 
staggered magnetization in both the cases, 
but with weaker effects for Ni than Zn. This 
is qualitatively similar to the present study 
on a much simpler system, i.e. in an ID 
chain, which confirms the validity of the ap- 
proach which consists in using magnetic or 
non-magnetic substitutions in the 2D system 
to probe the intrinsic magnetic correlations. 
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Figure captions 

FIG. 1: (a) Representation of the unit cell 
of Y2BaNiO"5 : the 89 Y nucleus is hyperfine- 
coupled to the Ni electrons in its two next 
nearest neighbor (n.n. n) chains through oxy- 
gens while its transferred hyperfine coupling 
to the Ni on the nearest neighbor (n.n.) chain 
is expected to be small, (b) When an impu- 
rity is substituted in one of the two n.n. n 
chains, 89 Y nuclei which couple to that chain 
are labelled according to their distance to the 
impurity as indicated in the figure. In Fig. 1 

(a) , only some of the oxygen atoms are rep- 
resented, for clarity. In Fig. 1 (b), only the 
Ni and Y sites are displayed. 

FIG. 2: Magnetic susceptibility versus 
temperature T for (a) Y2BaNii„ x Zn a; 05 and 

(b) Y 2 BaNii_ x Cu x 05 (0 < x < 0.02). The 
solid line follows Eq. 5 and the fitting was 
done using the procedure discussed in Sec. 
IVC, for the "pure" YBNO sample for T < 
150 K. In the insets, the representation vs 
1/T evidences the Curie law behavior at 
low temperatures. Comparison between non- 
annealed and annealed sample in the inset of 
(a) shows that annealing leads to a reduc- 
tion of the Curie term. In contrast, Zn or Cu 
substitution leads to an increase of the Curie 
term. 

FIG. 3: Comparison of the 89 Y NMR spec- 
tra of the non-annealed and annealed YBNO 
at 150 K (in the inset, the full line and in 
the main frame, a vertical zoom). The inten- 
sity of the satellite peaks reduce after vacuum 
annealing showing that the corresponding de- 
fects are eliminated by annealing. 

FIG. 4: 89 Y Knight shift of the main- 
line 89 -K" as a function of temperature T, for 
pure (present results and from Shimizu et 
al. (Ref. [14])) and substituted samples. 
S9 K(T), which measure the uniform suscep- 
tibility, remains unaffected by substitutions. 
The solid line represents a fit to Eq. 1 for the 
Cu doped sample and is discussed in detail 
in Sec. IVA. 

FIG. 5 (a) Expanded view of the satel- 
lites of the 89 Y spectra of Y 2 BaNii_, E Zn r 05 



(0 < x < 0.02) at 120 K. Intensity of these 
satellites increases with increasing Zn con- 
tent. The positions of the satellite peaks 
are independent of the Zn concentration, (b) 
Comparison of the satellites for 0.5% and 1% 
Cu substitution at the Ni site at 300 K. The 
satellite intensity increases linearly with Cu 
content. 

FIG. 6: 1 st and 3 rd neighbor satellite peaks 
(labelled T and '3') for Y2BaNio.9sZno.02O5 
and Y2BaNio.g8Cuo.02O5 at 200 K (solid lines) 
and fit of satellite 1 using an asymmetric 
Lorenzian lineshape (dotted lines). 

FIG. 7: Comparison of 89 Y NMR 
spectra of Y2BaNio.98Zno.02O5 and 
Y2BaNio.98Cuo.02O5 measured at T = 200 
K. The mainline position is the same for 
both substitutions. The satellites for the 
Cu substitution are less shifted than those 
for Zn. Inset: Position of the satellite 
peaks on the lower frequency side of the 
mainline for non-annealed pure YBNO, 
Y2BaNio.9sZno.02O5 and Y2BaNio.9sCuo.02O5 
at 200 K. 

FIG. 8: 89 Y NMR spectra of 
(a) Y2BaNio.9sZno.02O5 and (b) 
Y 2 BaNi .995Cuo.oo505 at different tem- 
peratures. Satellite peaks are indexed in 
accordance with their decreasing shift from 
the mainline and with the position of the 89 Y 
labelled in Fig. 1 (b). The peak originating 
from the Y neighbor to the Zn impurity (as 
discussed in Sec. IIIB) is marked by in (a). 

FIG. 9: Temperature dependence of 5ui, 
the frequency shift of the different satel- 
lites with respect to the central line for (a) 
YBNO:Zn and (b) YBNO:Cu. The open cir- 
cles in (a) represent the shifts obtained in 
Ref. [17] for YBNO:Mg. At lower temper- 
atures, the broadening of the NMR spectra 
restricted us from obtaining satellite shifts 
with high resolution. Indices represent the i th 
neighbor to the impurity in accordance with 
Fig. 1 (b). 

FIG. 10: 89 Y NMR shift of the central line 
K as a function of Xu computed by QMC us- 
ing different J values (with T as an implicit 
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parameter). The solid line represents the best 
linear fit and allows to determine both J and 
the hyperfme coupling constant Ahf (using 
the constraint K(0) = 370 ± 10 ppm as de- 
tailed in the text). 

FIG. 11: Site dependence of the impu- 
rity induced magnetic shift \5v>i\ at differ- 
ent temperatures for (a) YBNO:Zn and (b) 
YBNO:Cu. The log-scale representation ev- 
idences that \5vi\ is exponentially decaying. 
The solid lines are fits to the exponential re- 
lation represented by Eq. 2. 

FIG. 12: Comparison of £i mp (T) obtained 
from our experiments from Eq. 2 and Fig. 11 
for YBNO:Zn, YBNO:Cu , for YBNO:Mg by 
Tedoldi et al. (Ref. [17]), and for the infinite 
chain correlation length £(T) calculated by 
QMC by Kim et al. (Ref. [33]). 

FIG. 13: (a) Comparison at T = 100 K 
of the spatial dependence of the staggered 
magnetization 5 < S l z > / < S z > u for 
YBNO:Zn with the corresponding QMC sim- 
ulations for various J values, which shows 
a good agreement for J = 270 ± 10 K. (b) 
Comparison of the T variation of the magne- 
tization < S l z > of the 1 st and 2 nd neighbor 
(i — 1 and 2, respectively) with the QMC 
calculations. Filled symbols represent our ex- 
perimental data while open triangles and cir- 
cles represent our QMC results for J = 260 
K. The open stars are the results of classical 
Monte Carlo calculations (Ref. [19]). 

FIG. 14: Spatial dependence of the exper- 
imental ratio < 5S l z >c u / < SS l z >znof the 



staggered magnetization near Cu to that near 
Zn . The data are compared to QMC simu- 
lations for various values of the Heisenberg 
coupling J' between Cu and its n.n. neigh- 
bors Ni at T = 100 K. 

FIG. 15: The impurity induced total sus- 
ceptibility xfJ^ R {T) estimated from NMR 
data are plotted versus 1/T for YBNO:Zn!% 
(open triangles) and YBNO:Cui% (open 
squares). The solid lines represents the 
Curie-like fits to X^ R ( T )- Th e solid sym- 
bols represent the values obtained at 100 K 
from the corresponding QMC simulations. 

FIG. 16: Macroscopic susceptibility of 
pure YBNO and its decomposition into a sum 
of three terms : the QMC-computed (Hal- 
dane gap related) decrease of the uniform 
susceptibility for J = 260 K, a Curie term 
extracted from a fit to the data for T < 30 
K , and the residual part which cannot be 
accounted for by the two preceding forms. 

FIG. 17 : The residual term deduced from 
the macroscopic susceptibility analysis de- 
tailed in the text and in Fig. 16 are compared 
here for various samples. 

FIG. 18 : Variation of the Curie constant 
(jSquid ex t rac ted from the low-T susceptibil- 
ity measurements versus impurity concentra- 
tion for YBNO:Zn and YBNO:Cu. The linear 
fits (solid lines) were done without taking into 
account the Curie-term of the "pure" sample, 
as the use of substituent might reduce the 
concentration of native defects in the latter. 
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